%bandpass filter

function [filtered] = bpfil(x, low, high)

len = length(x)/2;
[b, a] = butter(1, [ (low/len), (high/len) ]);
filtered = filter(b, a, x);